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Abstract. We first recall the concept of Dyadosphere (electron-positron-photon plasma around a formed black holes) and 
its motivation, and recall on (i) the Dirac process: annihilation of electron-positron pairs to photons; (ii) the Breit- Wheeler 
process: production of electron-positron pairs by photons with the energy larger than electron-positron mass threshold; the 
Sauter-Euler-Heisenberg effective Lagrangian and rate for the process of electron-positron production in a constant electric 
field. We present a general formula for the pair-production rate in the semi-classical treatment of quantum mechanical 
tunneling. We also present in the Quantum Electro-Dynamics framework, the calculations of the Schwinger rate and effective 
Lagrangian for constant electromagnetic fields. We give a review on the electron-positron plasma oscillation in constant 
electric fields, and its interaction with photons leading to energy and number equipartition of photons, electrons and positrons. 
The possibility of creating an overcritical field in astrophysical condition is pointed out. We present the discussions and 
calculations on (i) energy extraction from gravitational collapse; (ii) the formation of Dyadosphere in gravitational collapsing 
process, and (iii) its hydrodynamical expansion in Reissner Nordstrom geometry. We calculate the spectrum and flux of photon 
radiation at the point of transparency, and make predictions for short Gamma-Ray Bursts. 
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INTRODUCTION 

Motivations and Dyadosphere. It is an one of most im- 
portant issues in modern physics to understand how grav- 
itational energy transforms to electromagnetic and rota- 
tional energies to during the process of gravitational col- 
lapses to black holes, in connection with observations. 
The primal steps toward the understanding of this issue 
are studies of electromagnetic properties of spinning and 
non-spinning black holes: (i) reversible and irreversible 
transformations - the Christodoulou-Ruffini formula, 
(ii) electron-positron pair-production in Kerr-Newmann 
geometry - the Damour-Ruffini proposal for Gamma 
Ray Bursts (GRBs), (iii) formation of electron-positron- 
photon plasma - Preparata-Ruffini-Xue Dyadopshere. 
Pair-production. The annihilation of electron-positron 
pair to two photons, and its inverse process - the pro- 
duction of electron-positron pair by the collision of 
two photons, as well as the electron-positron pair pro- 
duction from the vacuum in constant electromagnetic 
fields, were studied in quantum mechanics by Dirac, 
Breit, Wheeler, Sauter, Euler, Heisenberg respectively 
in 1930's, and Schwinger in Quantum Electro-Dynamics 
(QED)inl951. 

Nonuniform fields. It has been a difficult task to ob- 
tain the rate of electron-positron pair production in vary- 
ing electromagnetic fields in space and time. This issue 
has attracted much attention not only for its theoretical 
viewpoint, but also its possible applications in heavy-ion 
collisions and high-energy laser beams, as well as astro- 
physics. 



Plasma oscillation. A naive expectation is that such ex- 
ternal electric field rapidly vanishes for its source neu- 
tralized by electrons and positrons produced. However, 
the back reaction (screening effect) of electron-positron 
pairs on the external fields leads to the plasma oscilla- 
tion phenomenon: electrons and positron oscillating back 
and forth in phase with alternating electric field. Beside 
electron-positron pairs oscillating together with the elec- 
tric field, they interact with photons via the Dirac and 
Breit- Wheeler processes, and approach to a thermal con- 
figuration. 

Critical fields on the surface of massive nuclear cores. 
In ground based laboratories, it is rather difficult to built 
up electromagnetic fields at the order of the critical field 
value E c in macroscopic space-time scales. However, in 
the arena of astrophysics, supercritical electric fields are 
energetic-favorably developed on the surface of neutron 
star cores, due to strong, electroweak and gravitational 
interactions of degenerate nucleons and electrons. 
Dyadosphere formed in gravitational collapse. Initiat- 
ing with supercritical electric fields on the surface, gravi- 
tational collapses of nuclear massive cores and processes 
of pair production, annihilation and oscillation lead to the 
formation of high energetic and dense plasma of elec- 
trons, positrons and photons, Dyadosphere that we pro- 
posed in 1998. 

Hydrodynamic expansion after gravitational collapse. 
The adiabatic and hydrodynamic expansion of the 
electron-positron-photon plasma after gravitational 
collapse, up to the transparency to photons, account 
for daily observing phenomena of Gamma Ray Bursts 



(GRBs). 

Predications in connection with short Gamma ray 
Bursts. Armed with a complete knowledge of all these 
fundamental processes, we present our understanding 
on the genuine origin of GRB-phenomenon, and make 
some predictions in connection with observations of 
short GRBs. 



BASIC MOTIVATIONS AND 
DYADOSPHERE 

Energetics of Electromagnetic Black Holes. 

The process of gravitational collapse of a massive 
core generally leads to a black hole characterized by 
all the three fundamental parameters: the mass-energy 
M, the angular momentum L, and charge Q QJJ] . The 
phenomenon of gravitational collapse is crucial for the 
evolution of the system. Nonetheless in order not to 
involve its complex dynamics at this stage, we assume 
that the collapse has already occurred. Correspondingly 
a generally charged and rotating black hole has been 
formed whose curved space-time is described by the 
stationary Kerr-Newmann geometry in Boyer-Lindquist 
coordinates (t,r,6,<j>) 
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where A = r 2 



r 2 + a 2 cos 2 9, 



2Mr + a 2 + Q 2 and E 
a = L/M being the angular momentum per unit mass 
of the black hole. The Reissner-Nordstrom and Kerr 
geometries are particular cases for non-rotating a = 0, 
and uncharged 2 = 0, black holes respectively. The total 
energy in terms of the Coulomb and rotational energies 
is described by the Christodoulou-Ruffini mass formula 
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where M lv is the irreducible mass. The reversible (irre- 
ducible) process of the black hole, characterized by con- 
stant (increasing) irreducible mass, can (cannot) be in- 
verted bringing the black hole to its original state. En- 
ergy can be extracted approaching arbitrarily close to 
reversible processes which are the most efficient ones. 
Namely, from Eq. © it follows that up to 29% of the 
mass-energy of an extreme Kerr black hole (M 2 = a 2 ) 
stored in its rotational energy term 2 [f M . , whereas up to 



50% of the mass energy of an extreme EMBH with (Q — 

c 2q2 

M) stored in the electromagnetic energy term , can 
be in principle extracted. 

Vacuum polarization around an Electromagnetic Black 
Hole. 

It was pointed [3] that via Sauter, Heisenberg, Euler 
and Schwinger process, the electron-positron pair pro- 
duction occurring around an superritical Electromagnetic 
Black Hole (EMBH) is actually a very efficient almost 
reversible process of energy extraction, and extractable 
energy is up to 10 54 ergs that accounts for very energetic 
phenomenon of GRBs. In order to study the pair pro- 
duction in the Kerr-Newmann geometry, at each event 
(t,r,9,(j)) a local Lorentz frame is introduced, associated 
with a stationary observer & at the event (t,r,9,<p). A 
convenient frame is defined by the following orthogonal 
tetrad 
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In the so fixed Lorentz frame, the electric potential Ao, 
the electric field E and the magnetic field B are given by 
the following formulas, 

A = o)i 0) A", E^ffljjV* ^ = l -cofe a ^F yS . 

(7) 

One then obtains Aq = — Qr(LA) l < 2 , while the electro- 
magnetic fields E and B are parallel to the direction of 
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respectively. The spatial variation scale GM/c 2 of these 
background fields is much larger than the Compton 
wavelength h/m e c of the quantum field, then, for what 
concern pair production, it is possible to consider the 
electric and magnetic fields defined by Eqs. ( 18191 ) as con- 
stants in a neighborhood of a few wavelengths around 
any events (r,9,(p,t). Based on the equivalence princi- 
ple, the rate of pair-production process in a constant field 
over a flat space-time can be locally applied to the case 
of the curved Kerr-Newmann geometry: 
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where the critical field E c = m 2 c 3 /eh. It was assumed 
that electron and positron produced fly apart from each 
other, one goes inward to neutralize EMBHs and another 
goes to infinity. This view was fundamentally modified in 
Refs. |4j,|5|,|6j] by the novel concept of the Dyadosphere. 
Dyadosphere: electron-positron-photon plasma. 

We start with the Reissner-Nordstrtim black holes and 
consider a spherical shell of proper thickness 8 = ^ <C 
^ centered on the EMBH, the electric field is approx- 
imately constant in it. We can then at each value of the 
radius r model the electric field as created by a capacitor 
of width 8 and surface charge density o(r) = and 
express Eq. ( fTOt as, 
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where electric field E = Ana, O c = j^E c is the critical 
surface charge density. The pair creation process in these 
shells will continue until a value of the surface charge 
density reaches the critical value o c , and it takes 
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This time is so short that the light travel time is smaller 
or approximately equal to the width 8. Under these cir- 
cumstances the correlation between shells can be approx- 
imately neglected, thus we can justify the approximation 
of describing the pair creation process shell by shell. 

The Dyadosphere is composed by these shells from 
the horizon r + to r c / s , which is given by E(r c i s ) = E c and 
can be expressed as, 
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using the Planck charge q ? = {he) 1 1 2 and the Planck 
mass m ? = (hc/G) 1 ^ 2 , which clearly shows the hybrid 
gravitational and quantum nature of this quantity. The 
total number of shells is about (r^ — r+)/ and the total 
number of pairs, 
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We calculate the number and energy densities of pairs in 
the Dyadosphere 
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as shown in Figs.Q] and total energy is then 
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Due to the very large pair density given by Eq. (fT~5T > 
and to the sizes of the cross-sections for the process 
e + e~ ^7+7, the system is expected to thermalize to 
a plasma configuration for which 
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In Fig. [2] the total energy ( fT6T > and the average energy 

E to l _ 

per pair N e + e are shown in terms of EMBH's mass 

ji and charge Recently, the Dyadotorus: the plasma 
of electron-positron-photon created in Kerr-Newmann 
black holes is studied 07[]. 



PRODUCTION OF 
ELECTRON-POSITRON PAIRS 

Early quantum electrodynamics. 

We recall three results, which played a crucial role 
in the development of the Quantum Electro-Dynamics 
(QED). The first is the Dirac process of an electron- 
positron pair annihilation into two photons, 



e^+e ^yi + 7>, 
and the cross-section in the rest frame of electron: 
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where y= ^" + /m e c 2 ^> 1 is the energy of the positron 
and a = e 2 j [Antic) is the fine structure constant. The 
second is the Breit- Wheeler process of electron-positron 
pair production by two photons collision, 
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which is the inverse Dirac -process ( fT8l and the cross- 
section is related to ( fT8l by the ^fT-theorem , 
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FIGURE 1. The number-density n e + e - (r) (left) and average energy per pair in MeV (right) are plotted as a function of the radial 
coordinate for fi = M/Mq = 10 and % = Q/M = 1 (upper curve) and <^ = 0. 1 (lower curve). 




FIGURE 2. Left: Total energy of Dyadosphere as a function of EMBHs' mass and charge parameters n, Right: The average 
energy per pair is shown here as a function of the EMBH mass in solar mass units for t, = 1 (solid line), t, = 0.5 (dashed line) and 
<| = 0.1 (dashed and dotted line). 



where /3 is the relative velocity of electron and positron. 
The third is the vacuum polarization in external uniform 
electromagnetic field, studied by Heisenberg and Euler, 
following Sauter's work on quantum tunneling probabil- 
ity from negative energy states [see Fig. (fj)], 
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E c ~ 1.3 • 10 16 V/cm, 10 6 larger than the value required 
to ionize a hydrogen atom. Heisenberg and Euler ob- 
tained nonlinear Lagrangian from the Dirac theory, 
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and its series expansion in powers of a, 
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They found facts that A«5f e ff is a complex function of E 
and B, the imaginary part is associated with pair produc- 
tion when the electric field E > E c , and the vacuum be- 
haves as a dielectric and permeable medium in which, 
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where complex and jU,j are the field-dependent dielec- 
tric and permeability tensors of the vacuum. 
Quantum electrodynamics. 

The QED-Lagrangian describing the interacting sys- 
tem of photons, electrons, and positrons reads 
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where Jz?q e ^ is for free electrons, positron and photons. 
An external field A^ is incorporated by adding to the 
quantum field A^ in 
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The amplitude for the vacuum to vacuum transition in the 
presence ofA e : 
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The effective action as a functional of A e is: 

A,6/ eff [A e ] = -/ln(0|0). 
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Under the assumption that A e (x) varies smoothly over a 
finite spacetime region, we may define an approximately 
local effective Lagrangian A_S? e ff [A e (x)], 

Aj4 ff [A e ]~ [ d 4 xA5f ef[ [A e (x)} « VAfA.Sf eff [A e ], (30) 



where V is the spatial volume and time interval Af , over 
which the field is nonzero. The amplitude of the vacuum 
to vacuum transition (1281 has the form, 
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where vacuum-energy difference ASq — So(A e ) — So(0), 
and r is the vacuum decay rate. The probability that the 
vacuum remains as it is in the presence of the external 
field A e is 

|(0|0)| 2 = e- 2ImA ^[ Ae ]. (32) 

This determines the decay rate of the vacuum caused by 
the production of electron and positron pairs: 
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We calculate the imaginary part d33l > and reproduce the 
Schwinger formula, 
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where e 2 — j3 2 = E 2 — B 2 and ej3 = EB. In addition, we 
calculate the real part d34l 
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and Ei(z) is the exponential-integral function. In the 
weak-field expansion, we obtain Eq. (124-b . In the strong- 
field expansion, we obtain, 
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FIGURE 3. In presence of a strong enough electric field the boundaries of the classically allowed states ("positive" or "negative") 
can be so tilted that a "negative" is at the same level as a "positive" (level crossing). Therefore a "negative" wave-packet from the 
left will be partially transmitted, after an exponential damping due to the tunneling through the classically forbidden states, as s 
"positive" wave-packet outgoing to the right, /i is particle's mass, eV potential energy and to energy. 



In the case E 3> 1, B = and m — 0, we obtain, 
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Nonuniform electric fields. 

Let the field vector E(z) point in the z-direction. The 
one-dimensional electric potential Aq(z) = — J z dz'E{z') 
and the positive and negative continuum energy-spectra 
are 



S ± =±^(cp z ) 2 + c 2 pl + (m e c 2 ) 2 + V(z), (41) 

where p z is the classical momentum in z-direction, px 
transverse momenta, and V (z) — eAo(z) potential energy. 
The crossing energy-levels between two energy-spectra 
S± (f4Tb appear, e = $+ = £-. The probability amplitude 
for quantum tunneling process can be estimated by a 
semi-classical calculation using WKB method (see e.g. 
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Changing the variable of integration from z to y(z), 
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where E(y) = E[z(y)]/Eo and z(y) is the inverse function 
of Eq. (l44l . The flux density of virtual particles attempt 
for tunneling at z- is 
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and the energy-variation t/e = |eE(z_)|dz. Using 
Eqs. ( 145|46t and expanding up to /(m 2 c 2 ), we obtain 
the WKB-rate of pair production per unit volume at 
given crossing energy-level e(z-) is 
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where V = V±dz,V± = J dxdy; D s = 2 for a spin- 1/2 
particle and D s = 1 for spin-0. The G and g functions 
are 
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FIGURE 4. Positive and negative energy-spectra £±(z) of 
Eq. J4U in units of m e c 2 , with p z = = as a function of 
x = z/l for the Sauter potential V± (z) J51t for <j = 5. 



Eq. ( f4Tb gives the semi-classical WKB-rate of pair- 
production per unit volume for any one-dimensional 
electric field E(z) and potential V(z) [10], provided 
crossing energy-levels e between negative and positive 
energy-spectrum occur. 

We apply our formula ( 147148b to a uniform field case, 
obtain G = l,g = G/2, 
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which is independent of crossing energy-levels e and 
coordinate z. gives the Sauter factor (|22| | and Heisenberg- 
Euler pref actor obtained from d23t . 
Sauter electric field. 

As an example, we consider the nonuniform Sauter 
electric field localized within finite slice of space of 
width t in the xy-plane 11011 . Electric field E — E(z)z and 
potential energy V(z) are given by 

E(z) = £ /cosh 2 (z/£) , (50) 
V(z) = am e c 2 twh(z/t), (51) 



where 
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and see Fig.|4] Using our formula ( 147148b and integrating 
over e(z-), we approximately obtain 
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The comparison between pair-production rates in un- 
bound uniform and bound nonuniform fields is given by 
the ratio /? rate of Eq. d49l and Eq. d53l 
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In Fig. ©, it is shown that the pair-production rate in 
the Sauter field becomes smaller as the confining size 



of the field becomes smaller. In Ref. ifioll . we present 
detailed calculations and discussions of pair-productions 
rate in various cases of nonuniform electric fields: the 
Sauter and Coulomb fields, as well as fields E(z) ~ z and 
E(z) t^0,z>0. 



PLASMA OSCILLATIONS OF 
ELECTRON-POSITRON PAIRS IN 
ELECTRIC FIELDS 

We have discussed the Sauter-Heisenberg-Euler- 
Schwinger process for electron-positron pair production. 
However, we neglect the following dynamics: 

1 . the back reaction of pair production on the external 
electric field; 

2. the screening effect of pairs on the external electric 
field strengths; 

3. the motion of pairs and their interactions. 

When these dynamics are considered, a phenomenon 
of electron-positron oscillation, plasma oscillation, takes 
place. We quantitatively discuss this phenomenon by 
using the relativistic Boltzmann-Vlasov equations [11] 

dtfe + eE ■ V p / e = (E,p) - % (f ,p) , (55) 
a/ r = 2Sf r (f,k), (56) 

where f e (fy) is spatially independent distribution func- 
tion of electrons (photons) in phase space; and the ho- 
mogeneous Maxwell equations, 



^E = -j„ (E)-J c (t), 



(57) 



where j p is the polarization current and j,- conduction 
current. The terms ^ (f,p) and Ify (f,k) stand for colli- 
sions between electrons, positrons and photons. 5? (E,p) 
is related to the pair-production rate d49l . 
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where [l±2/ e ] accounts for the Bose enhancement(+) 
and Pauli blocking (-). 

These Equations 051157] ) are integrated with the fol- 
lowing initial conditions of |E| =Eq = 9E ( and null den- 
sities of electrons, positron and photons. The results of 
the numerical integration in units of the Compton time 
%c and length Ac are shown in Fig. [6] at early times, 

1 . the electric field does not abruptly reach the equi- 
librium value but rather oscillates with decreasing 
amplitude; 

2. electrons and positrons oscillates in the electric field 
direction, reaching ultra relativistic velocities; 




FIGURE 5. For E = E c , a = ijXc <HJ is the spatial size where electric field E ^ 0. The ratio ^ rate J54t is plotted as function 
of a in the left figure. The number of pairs created in Compton area and time as functions of a (up curve for Schwinger constant 
field l |49t and low one is nonuniform Sauter field l |53t ) in the right figure. 



3. the role of the e + e <=^ 77 scatterings is marginal in 
the early time of the evolution. 

At late times the system is expected to relax to a plasma 
configuration of thermal equipartition with the asymp- 
totic behavior: 

1. the electric field is screened to about the critical 
value: E ~ E c for t ~ 10 3 - 10 4 T C > T C ; 

2. the initial electromagnetic energy density is dis- 
tributed over electron-positron pairs and photons, 
indicating energy equipartition; 

3. photons and electron-positron pairs number densi- 
ties are asymptotically comparable, indicating num- 
ber equipartition. 



Note that we show 111211 that such phenomenon of plasma 
oscillation occurs not only for strong fields, but also for 
weak fields, in addition, a detailed study of thermaliza- 
tion of electrons-positrons-photons plasma is given in 
Ref. fl^l . The thermalized plasma starts hydrodynam- 
ical expansion described by hydrodynamic equations 

CUBE. 



SUPER CRITICAL FIELD ON THE 
SURFACE OF MASSIVE NUCLEAR 
CORES 

Electromagnetic properties of massive nuclear cores 

Based on numerical lfl9ll and analytical Ir20l I2H1 ap- 
proaches to the relativistic Thomas-Fermi theory, we 
study electron configurations and electromagnetic prop- 
erties of massive nuclear cores of mass number A and 
radius R, 



Mpi 



anck 



m N 



AO 51 , R c ~—N l J 3 ~10km,(59) 
m K c 



where Mp\ anc ^,mN,m K are Planck, nucleon and pion 
masses, with the global neutrality condition: the same 
proton and electron numbers N p = N e . We show that 



close to core surface, it exists supercritical electric field 
E > E c , prove that this configuration is stable and en- 
ergetically favorable against the configuration with the 
local neutrality condition: the same proton and electron 
densities n e (x) — n p (x), usually adopted. 

The Thomas-Fermi theory for the electrostatic equi- 
librium of electron distributions n e (r) around extended 
nuclear cores can be described as follow. Degenerate 
electron density n e (r), Fermi momentum Pf and Fermi- 
energy S e {Pf ) are related by 



F\3 



n e {r) = 



(P[) 
3n 2 h 



3 ' 



(60) 



S e {K) = [(P[c) 2 +m 2 e c 4 ] l / 2 -m e c^V conl (r), 

where V cou \(r) is Coulomb potential energy. The electro- 
static equilibrium of electron distributions is determined 
by 



(61) 



which means the balance of electron's kinetic and poten- 
tial energies in Eq. ( |60l > and degenerate electrons occupy 
energy-levels up to +m e c 2 . Eqs. ( 16016 It give: 



n e (r) 



[V 2 Jr)+2m e c 2 V coul (r)fl62) 



3n 2 (ch)' : 



The Gauss law leads the following Poisson equation and 
boundary conditions, 

AVcouiM = 4na [n p (r) - n e {r)] ; 
KouiH = 0, VW(0) = finite. (63) 

Degenerate proton and densities n pM ir) are constants 
inside core r < R c and vanishes outside the core r > R c . 



n p , n (r) = 



'"' nvJ 3;rV 

^ _ UuF A 2, m 2 411/2. _ „2 



(64) 



%AK) - l(Pp. n cr+m 2 M c«Y/ z -m p , n c z +V coal {r)8 p , 



1 




-10 I 1 ' ' 1 ' ' ' 1 1000 2000 3000 4000 5000 6000 7000 




FIGURE 6. In left figure: We plot for / < 150tc, from the top to the bottom panel: a) electromagnetic field strength; b) electrons 
energy density; c) electrons number density; d) photons energy density; e) photons number density as functions of time. The right 
figure: We plot for t < 7000tc as the same quantities as in left. 



where P* n i<Sp, n {Pe) are Fermi momenta, energies of 
protons and neutrons, and S p indicates V cou \(r) for pro- 
tons only. Neutrinos assumed to escape from massive 
cores, the energetic equation for the j3 —equilibrium of 
neutrons, protons and electrons is 

g n {f n )=£ p {Pl) + £ e {P F e \ (65) 

which gives the relationship between the neutron, proton 
and mass numbers N„,N P ,A = N„ + N p . We integrate 
these equations are integrated and show results in Fig. (jTJ 
The configuration is electrostatic stable, since the mean 
repulsive energy is much smaller than mean gravitational 
binding 0.1 M^c 2 for protons in the surface layer. 

GEOMETRY OF GRAVITATIONALLY 
COLLAPSING CORES 

The Tolman-Oppenheimer-Snyder solution. 

Oppenheimer and Snyder first found a solution of the 
Einstein equations describing the gravitational collapse 
of spherically symmetric star of mass greater than ~ 
O.7M . In this section we briefly review their pioneering 
work as presented in Ref. [22]. 

In a sperically symmetric space-time they can be 
found coordinates (t,r,9,(f>) such that the line element 



takes the form 

ds 2 = e v dt 2 - e X dr 2 - r 2 d& 2 , (66) 

dD 2 = dO 2 + sin 2 6d(j) 2 , v = v(t,r), X = X(t,r). How- 
ever the gravitational collapse problem is better solved in 
a system of coordinates (i,R, 9,(j)) which are comoving 
with the matter inside the star. In comoving coordinates 
the line element takes the form 

ds 2 = dx 2 - e a dR 2 - e m dQ. 2 , 

CO = W(z,R), (0 — G>(t,R). Einstein equations read 

8kTI = e - u> -e- c ^ + (0+l6) 2 (67) 
87rT 2 2 = 87rT 3 3 = (2»" + co' 2 - a' co') 

+ ±(2d + 6 2 + 2cb + co 2 + da>) (68) 

8ttT 4 4 = e- a - e-° (ft)" + f ft)' 2 - ^f) + % + f 

(69) 

87te a T^ = -8xT?= ift)'(ft>-(7) + aV. (70) 

Where T^y is the energy-momentum tensor of the stel- 
lar matter, a dot denotes a derivative with respect to T 
and a prime denotes a derivative with respect to R. Op- 
penheimer and Snyder were only able to integrate Eqs. 
(I67li-(l70"ii in the case when the pressure p of the stellar 
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FIGURE 7. Potential energy — V C oul( r ) (left) an d electric field E(r) ~ — V c ' oul (r) (right) are plotted as a function of (r- 

w,)/(ioa.-:. 



matter vanishes and no energy is radiated outwards. In 
the following we thus p = 0. In this hypothesis 



tv 1 = r 2 2 = r 3 3 = t} = T? = 0, 



where p is the comoving density of the star. Eq. ( |70| ) was 
first integrated by Tolman in Ref . J23ll . The solution is 



e° =e m CO ,2 /4f{R), 



(71) 



where / = f(R) is an arbitrary function. In Ref. l22tl was 
studied the case f(R) = 1. The hypothesis f(R) = 1 will 
be relaxed in the case of a shell of dust. Using Eq. ( TTTT > 
into Eq. d67l ) with f(R) = 1 gives 



0, 



which can be integrated to 

e a = {Fx + Gfl\ 



(72) 



(73) 



where F — F(R) and G = G(R) are arbitrary function. 
Using Eq. (|7TT i into Eq. d68l > gives Eq. (l72l again. From 
Eqs. (l69l >. (l7ll and (l73l the density p can be found as 



-l 



87Tp = |(T+f) '(T+^j . (74) 

There is still the gauge freedom of choosing R so to have 

G = R'I\ 

Moreover, it can be freely chosen the initial density 
profile, i.e., the density at the initial time T = 0, po = 
Po(R). Eq. d74l i then becomes 

FF' = 9nR 2 po(R) 

whose solution contains only one arbitrary integration 
constant. It is thus seen the the choice of Oppenheimer 
and Snyder of putting f(R) = 1 allows one to assign only 
a one-parameter family of functions for the initial values 
Po = Po{R) of p. However in general one should be able 



to assign the initial values of p arbitarrily. This will be 
done in the following section in the case of a shell of 
dust. 

Choosing, for instance, 



const > if R < Rb 
if R>Rb 



Rb being the comoving radius of the boundary of the star, 
gives 




where r+ = 2M is the Shwarzschild radius of the star. 

We are finally in the position of performing a co- 
ordinate transformation from the comoving coordinates 
(T,R,9,<p) to new coordinates (t,r,6,(j>) in which the 
line elements looks like d66l i. The requirement that the 
line element be the Schwarzschild one outside the star 
fix the form of such a coordinate transformation to be 



t = 




2 -l/2 f/? 3/2 

3 r + y K b 



. r ,y/2 ) _ 2r+y i/2. 



r 3 / 2 )-2(rr+) 1 / 2 + r H 



log 



,l/2 + 1/2 



r l/2-r 



where the first line for if R < Rb, the second line for 
if R>R b and 



R b r 
r+R- 



Gravitational collapse of charged and uncharged shells. 

It is well known that the role of exact solutions has 
been fundamental in the development of general relativ- 
ity. In this section, we present here exact solutions for 
a charged shell of matter collapsing into an Electromag- 
netic Black Hole (EMBH). Such solutions were found 



in Ref. 114711 and are new with respect to the Tolman- 
Oppenheimer-Snyder class. For simplicity we consider 
the case of zero angular momentum and spherical sym- 
metry. This problem is relevant for its own sake as an 
addition to the existing family of interesting exact solu- 
tions and also represents some progress in understanding 
the role of the formation of the horizon and of the irre- 
ducible mass discussed in Ref. II4611 . It is also essential 
in improving the treatment of the vacuum polarization 
processes occurring during the formation of an EMBH 
discussed in Ref. 1 1311 . As we already mentioned, both of 
these issues are becoming relevant to explaining gamma 
ray bursts, see e.g. l24T,l55ll25ll26| | and references therein. 

W. Israel and V. de La CruzlH |H] showed that the 
problem of a collapsing charged shell can be reduced 
to a set of ordinary differential equations. We recon- 
sider here the following relativistic system: a spherical 
shell of electrically charged dust which is moving radi- 
ally in the Reissner-Nordstrom background of an already 
formed nonrotating EMBH of mass M\ and charge Q\, 
with Q\ < Mi. The Einstein-Maxwell equations with a 
charged spherical dust as source are 



V [M F vp] 



r (d) 



■ 4nf 



where 



rH 
';UV 

T (d) 



35t (F/Fpv - 



pa) 



SltflUvi J 



aur 



(75) 



(76) 



Here Tj^j , 7^y" lj and jf are respectively the energy- 
momentum tensor of the dust, the energy-momentum 
tensor of the electromagnetic field F^ v and the charge 
4— current. The mass and charge density in the comoving 
frame are given by e, a and u a is the 4-velocity of the 
dust. In spherical-polar coordinates the line element is 



ds 2 = g fiv dx' i dx v = -e v{r ^dt 2 +e X{a Ur 2 

(77) 

where dQ 2 = d6 2 + sin 2 9d(j) 2 . 

We describe the shell by using the 4-dimensional 
Dirac distribution 5' 4 ' normalized as 



J S< 4 > (*,*') y/=g?X = 1 



(78) 



where g — det ||^v|[- We then have 

£ (x) = Mo j <5 (4) (x,x ) r 2 dxdQ., (79) 
o (x) = Q f <5 (4) (x,x ) r 2 dxdQ.. (80) 



Mq and Qo respectively are the rest mass and the charge 
of the shell and T is the proper time along the world 



surface S : xo — xo(t,£1) of the shell. S divides the 
space-time into two regions: an internal one and an 
external one As we will see in the next section for 
the description of the collapse we can choose either 
or ^# + . The two descriptions, clearly equivalent, will be 
relevant for the physical interpretation of the solutions. 
Introducing the orthonormal tetrad 



(0) 



ft) 



(3) 



fl /2 dt, 

rsin 9d(j), 



CO 



= f ± l/2 dr 7 CO^^rdO, 
(81) 



we obtain the tetrad components of the electric field 



S = £ ft) (1) = 



4 (B\) outside the shell 
%■ oL 1 ' inside the shell 



(82) 



where Q = Qo + Qi is the total charge of the system. 
From the G tt Einstein equation we get 



ds 



-f+dt 2 +f+ l dr 2 + r 2 dQ. 2 



-f-dti+fI l dr z + r z dQ. z 



outside the shell 
inside the shell 
(83) 

where /+ = 1 - ^ + /_ = 1 - 2£L + ^ and t- and 
t + are the Schwarzschild-like time coordinates in 
and respectively. Here M is the total mass-energy of 
the system formed by the shell and the EMBH, measured 
by an observer at rest at infinity. 

Indicating by ro the Schwarzschild-like radial coordi- 
nate of the shell and by to± its time coordinate, from the 
G tr Einstein equation we have 



Mo 

2 



(84) 

The remaining Einstein equations are identically satis- 
fied. From d84l l and the normalization condition u^u^ = 
— 1 we find 



(85) 

% = i^f M _ MlT ^_ ei_eieoy (86) 

(It M Q f±{r ) V 1 T 2i'fl 2r r J v ' 



We now define, as usual, r±= M ± \J M 2 — Q 2 : when 
Q < M, r± are real and they correspond to the horizons of 
the new black hole formed by the gravitational collapse 
of the shell. We similarly introduce the horizons — 

Mi ± ^M\-Q\ of the already formed EMBH. From 
we have that the inequality 



61 go 



>0 



(87) 



holds for r > r + if Q < M and for r > r\ if Q > M 
since in these cases the left hand side of d84l ) is clearly 
positive. Eqs. ([85]) and ([86]) (together with ([83]), d82b ) 
completely describe a 5-parameter (M, Q, M\, Q\, Mq) 
family of solutions of the Einstein-Maxwell equations. 

For astrophysical applications 1 13] the trajectory of the 
shell ro = ro (fo+) is obtained as a function of the time 
coordinate fo+ relative to the space-time region ^# + . In 
the following we drop the + index from fo+- From 
and (f86l> we have 



di'O — dr dl u£,/o2 _ F 

dt ~ dx dt ~ zn n v ii r ' 



(88) 



where 



F=f + (n>) = l-f 



ii = r 



2M r 



M-M l 
M 



(89) 



Since we are interested in an imploding shell, only the 
minus sign case in d88l will be studied. We can give the 
following physical interpretation of r. If M — M\ > Mq, 
r coincides with the Lorentz y factor of the imploding 
shell at infinity; from d88l it satisfies 



fdr 
\<"o 



> 1. 



(90) 



When M — M\ < Mo then there is a turning point rg, 
defined by ^ =0. In this case F coincides with 
the "effective potential" at r ( j : 



r=V/-(^)+V(-§ + l 



QiQo 



The solution of the differential equation 



by: 



dto = — 



=dr - 



<1. (91) 
is given 

(92) 



The functional form of the integral d92"] > cru- 
cially depends on the degree of the polynomial 
P( r o) = r o(^ 2 — which is generically two, but 
in special cases has lower values. We therefore distin- 
guish the following cases: 

1. M = Mq+M 1 ; Qi =M\\ Q = M: P(r ) is equal to 
0, we simply have 



ro(to) = const. 



(93) 



2. M = Mo+M x ;M 2 -Q 2 =M 2 -Q\;Q^M:P(ro) 



is a constant, we have 



const 



rt log ( 



2^/m 2 -q 2 



M 



3. M = M +M l ;M 2 -Q 2 ±M\- Q\: P (r ) is a first 
order polynomial and 



to = const + 2rov Q 2 — F 



Mm 



3(M 2 -Q 2 -M 2 +Q]) 
(M 2 +g 2 - Q\ ) 2 -9MM (Ml +Q 2 - Q\ ) + 1 2M 2 M 2 +2g 2 M 2 

3(M 2 -e 2 -M 2 +e 2 ) 2 

1 
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(95) 



where Q.± = Q. (r±). 
4. M Mo +M\\ P (ro) is a second order polynomial 
and 



to = const - 



■r 2 log 



r 2 _ log 



l 



\ 2YJM 2 -Q 2 - 

= [ r 2 -i r °vfl 2 -f 



2V^ 2 -e 2 

roVn 2- ? r 2 (n 2 -F)+r 2 n 2 + -(r 2 -i)( ro -r + ) 2 " 

D- r + 2(r -r + Wn 2 -F 

roVtf^F , r 2 (a 2 -F)+r 2 n 2 -(r 2 -l)(r -r-) 2 " 
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[2MM (2r 3 -3r)+M 2 +S 2 -e 2 ]V M2 -2 2 , 
2M ( r 2 - 1 ) ro - (M 2 + Q 2 - Q 2 ) T+2M M 



(96) 



(r + 2)r + r2log(2^±) 
(94) 



In the case of a shell falling in a fiat background {M\ = 
Q\ = 0) it is of particular interest to study the turning 
points rj5 of the shell trajectory. In this case equation d85b 
reduces to 

(S) 2 = 4(*+|-f) ! -i- m 

Case (2) has no counterpart in this new regime and 
Eq. (f87l > constrains the possible solutions to only the 
following cases: 

1. M = M ; Q = Mo- r = r (0) constantly. 

2. M = Mo; Q < Mo- There are no turning points, the 
shell starts at rest at infinity and collapses until a 
Reissner-Nordstrom blac k-hole is f ormed with hori- 
zons at rg = r± = M ± y/ M 2 — Q 2 and the singular- 
ity in ro = 0. 

3. M Mq. There is one turning point rg. 

(a) M < Mq, then necessarily is Q < Mo- Posi- 
tivity of rhs of d97l ) requires ro < rg, where 

* l O 2 — Mx ... . . „, 

r o = I m-Mq ls unique turning point. Then 
the shell starts from rg and collapses until the 
singularity at ro = is reached. 

(b) M > Mq- The shell has finite radial velocity at 
infinity. 

i. Q < Mq. The dynamics are qualitatively 
analogous to case (2). 



ii. Q > Mq. Positivity of the rhs of (J97j and 
requires that ro > r^, where r^ — 



3 m-Mo ■ ^ ne srie U starts from infinity and 
bounces at ro = rL reversing its motion. 

In this regime the analytic forms of the solutions are 
given by Eqs. d95l > and ( |%1 >, simply setting M\ = Q\ = 0. 

Of course, it is of particular interest for the is- 
sue of vacuum polarization the time varying electric 
field <=y = % on the external surface of the shell. 

In order to study the variability of S rQ with time 
it is useful to consider in the tridimensional space 
of parameters (Vo , f o , <? ro ) the parametric curve ^ : 

ro = A, to = fo(A), (o ro = j^j ■ In astrophysical ap- 
plications lfl3ll we are specially interested in the family 
of solutions such that 4? is when rn = °° which im- 
plies that T = 1. In Fig. [8] we plot the collapse curves 
in the plane (?o,ro) for different values of the parameter 
% = %, < % < I, The initial data are chosen so that 
the integration constant in Eq. d95l l is equal to 0. In all 
the cases we can follow the details of the approach to 
the horizon which is reached in an infinite Schwarzschild 
time coordinate. 
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FIGURE 8. Collapse curves in the plane (T,R) for M = 
2OM and for different values of the parameter t, . The asymp- 
totic behavior is the clear manifestation of general relativistic 
effects as the horizon of the EMBH is approached 

In Fig. [9] we plot the parametric curves 1? in the space 
(ro,?o, S'rg) for different values of % . Again we can follow 
the exact asymptotic behavior of the curves ^ , <o ro reach- 



ing the asymptotic value %■ . The detailed knowledge of 

this asymptotic behavior is of great relevance for the ob- 
servational properties of the EMBH formation, see e.g. 
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FIGURE 9. Electric field behaviour at the surface of the shell 
for M = 20Mq and for different values of the parameter %. 
The asymptotic behavior is the clear manifestation of general 
relativistic effects as the horizon of the EMBH is approached 



Irreducible mass of an EMBH and energy extraction 
processes 

The main objective of this section is to clarify the 
interpretation of the mass-energy formula [2] for an 
EMBH. For simplicity we study the case of a nonro- 
tating EMBH using the results presented in the previ- 
ous section. As we saw there, the collapse of a nonrotat- 
ing charged shell can be described by exact analytic so- 
lutions of the Einstein-Maxwell equations. Consider to 
two complementary regions in which the world surface 
of the shell divides the space-time: .M- and ^#+. They 
are static space-times; we denote their time-like Killing 
vectors by ^ and respectively. is foliated by the 
family {E, + : t + = t } of space-like hypersurfaces of con- 
stant t + . 

The splitting of the space-time into the regions jM- 
and allows two physically equivalent descriptions 
of the collapse and the use of one or the other depends 
on the question one is studying. The use of proves 
helpful for the identification of the physical constituents 



of the irreducible mass while jM+ is needed to describe 
the energy extraction process from EMBH. The equation 
of motion for the shell, Eq. [85] reduces in this case to 



(^) 2 =( 



in , 



and 



(<£) 2 =(M-g-ii) -Ml h 
y. The constraint l87lbecomes 



M 



Q 1 



>0. 



(98) 



(99) 



(100) 



Since is a flat space-time we can interpret — ^ in 

d98l as the gravitational binding energy of the system. ^ 
is its electromagnetic energy. Then Eqs. j98l , d99b dif- 
fer by the gravitational and electromagnetic self-energy 
terms from the corresponding equations of motion of a 
test particle. 

Introducing the total radial momentum P = Mqu 1 ' = 
Mq^j of the shell, we can express the kinetic energy of 
the shell as measured by static observers in as T = 

-M u^ -M = yJP* + M$-M . Then from Eq. 
we have 



2r 



St 
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pa. 



-M 2 = 



M Q - 



' 1 2r + 2r ' 



(101) 

where we choose the positive root solution due to the 
constraint ( 1 100b . Eq. (1101b is the mass formula of the 
shell, which depends on the time-dependent radial coor- 
dinate ro and kinetic energy T. If M > Q, an EMBH is 
formed and we have 



M = M + T + 



2r H 



el 

2r+ ' 



(102) 



where T+ = T (r+) and r+ = M + ^M 2 - Q 2 is the 
radius of the external horizon of that 



M=M„ 



St 



so it follows that 



Mir = M - 



(103) 



(104) 



namely thatMi r is the sum of only three contributions: the 
rest mass Mo, the gravitational potential energy and the 
kinetic energy of the rest mass evaluated at the horizon. 
M; r is independent of the electromagnetic energy, a fact 
noticed by Bekenstein J29ll . We have taken one further 
step here by identifying the independent physical con- 
tributions to Ma. This will have important consequences 
for the energetics of black hole formation (see 114611 ). 



Next we consider the physical interpretation of the 
ectromagnetic term v 
uating the Killing integral 



electromagnetic term which can be obtained by eval- 



j^T^diy = J" Sdr^d cos 9 
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(105) 



where is the space-like hypersurface in described 
by the equation t+ = t = const, with dT. v as its surface 
element vector. The quantity in Eq. (1105b differs from 
the purely electromagnetic energy 



1+ Jr 



(106) 



where — f + ' ^ is the unit normal to the integra- 
tion hypersurface and g rr = f + . This is similar to the 
analogous situation for the total energy of a static spher- 
ical star of energy density e within a radius ro, m {tq) = 
4-n fo° dr r 2 e, which differs from the pure matter energy 
m p( r o) = 47t /q° dr^/g^r^e by the gravitational energy 
(see i36Tl ). Therefore the term ^ in the mass formula 
dlOU is the total energy of the electromagnetic field and 
includes its own gravitational binding energy. This en- 
ergy is stored throughout the region ./#+, extending from 
ro to infinity. 

We now turn to the problem of extracting the electro- 
magnetic energy from an EMBH (see [2]). We can distin- 
guish between two conceptually physically different pro- 
cesses, depending on whether the electric field strength 
£ = % is smaller or greater than the critical value <§ c . 



The maximum value , 



_ 2 



of the electric field around 



an EMBH is reached at the horizon. We then have the 
following: 

1. For §+ < # c me leading energy extraction mecha- 
nism consists of a sequence of descrete elementary 
decay processes of a particle into two oppositely 
charged particles. The condition <S + < S z implies 



GM/c 



s/GM 



\VGme) 



-1 



1Q -6 M_ 



where the first line is for if M- < 10 6 , the second 
line for if M- > 10 6 and Ac is me Compton wave- 
length of the electron. Denardo and Ruffini i37ll 
and Denardo, Hively and Ruffini ll38ll have defined 
as the effective ergosphere the region around an 
EMBH where the energy extraction processes oc- 
cur. This region extends from the horizon r+ up to a 



radius 



The number of particles created [4] is then 
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(107) 

The energy extraction occurs in a finite number 
AfpD of such discrete elementary processes, each one 
corresponding to a decrease of the EMBH charge. 
We have 

(108) 



Since the total extracted energy is (see Eq 
ptot _ Q 1 

c — 2. 



we obtain for the mean energy per ac- 



celerated particle (£) PD — fj^ 



(F) - -21 

\ C /PD — 2r + 



which gives 
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(109) 
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77*e theorem of maximum energy extraction from 
gravitational collapse. 

One of the crucial aspects of the energy extraction 
process from an EMBH is its back reaction on the 
irreducible mass expressed in [2]. Although the en- 
ergy extraction processes can occur in the entire ef- 
fective ergosphere defined by Eq. (11071 1. only the 
limiting processes occurring on the horizon with 
zero kinetic energy can reach the maximum effi- 
ciency while approaching the condition of total re- 
versibility (see Fig. 2 in H for details). The farther 
from the horizon that a decay occurs, the more it 
increases the irreducible mass and loses efficiency. 
Only in the complete reversibility limit [2] can the 
energy extraction process from an extreme EMBH 
reach the upper value of 50% of the total EMBH 
energy. 

2. As we already discussed, for S + > S c the leading 
extraction process is the collective process based on 
the electron-positron plasma generated by the vac- 
uum polarization. The condition <? + > S c implies 



GM/c 2 I e \ 
hi \ \[Gm e J 



-1 



2-lQ- 6 ^- < 1. (Ill) 



This vacuum polarization process can occur only 
for an EMBH with mass smaller than 2 • 1O 6 M . 
The electron-positron pairs are now produced in the 
dyadosphere of the EMBH. We have 
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The total energy stored in the dyadosphere is Q] 
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which can be also rewritten as 
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We stress again that the vacuum polarization around 
an EMBH has been observed to reach the maximum 
efficiency limit of 50% of the total mass-energy of 
an extreme EMBH (see e.g. ]4[]). The conceptual 
justification of this result needs, however, the dy- 
namical analysis of the vacuum polarization process 
during the gravitational collapse and the implemen- 
tation of the screening of the e + e~ neutral plasma 
generated in this process. This analysis conceptu- 
ally validates the reversibility of the process and is 
given in the next chapter. 

Let us now compare and contrast these two processes. 
We have 
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Moreover we see (Eqs. dl 101 ). ( II 161 >) that (E) PD is in the 
range of energies of UHECR (see [30] and references 
therein), while for | ~ 0.1 and M ~ 1OM , (E) d is in 
the gamma ray range. In other words, the discrete par- 
ticle decay process involves a small number of particles 
with ultra high energies (~ 10 21 eV), while vacuum polar- 
ization involves a much larger number of particles with 
lower mean energies (~ lOMeV). 

The new conceptual understanding of the mass for- 
mula has important consequences for the energetics of 
a black hole. The expression for the irreducible mass in 
terms of its different physical constituents (Eq. ( 1104b ) 
leads to a reinterpretation of the energy extraction pro- 
cess during the formation of a black hole as expressed 



in 04611 . It will certainly be interesting to reach an under- 
standing of the new expression for the irreducible mass 
in terms of its thermodynamical analogues. 
The theorem of the maximum energy extraction from 
gravitational collapse. 

In this section we turn to the following aim: pointing 
out how formula 11041 for the M m leads to a deeper 
physical understanding of the role of the gravitational 
interaction in the maximum energy extraction process 
of an EMBH. This formula can also be of assistance 
in clarifying some long lasting epistemological issue 
on the role of general relativity, quantum theory and 
thermodynamics. 

It is well known that if a spherically symmetric mass 
distribution without any electromagnetic structure under- 
goes free gravitational collapse, its total mass-energy M 
is conserved according to the Birkhoff theorem: the in- 
crease in the kinetic energy of implosion is balanced 
by the increase in the gravitational energy of the sys- 
tem. If one considers the possibility that part of the ki- 
netic energy of implosion is extracted then the situation is 
very different: configurations of smaller mass-energy and 
greater density can be attained without violating Birkhoff 
theorem. 

We illustrate our considerations with two examples. 
Concerning the first example, it is well known from the 
work of Landau [31] that at the endpoint of thermonu- 
clear evolution, the gravitational collapse of a spherically 
symmetric star can be stopped by the Fermi pressure of 
the degenerate electron gas (white dwarf). A configura- 
tion of equilibrium can be found all the way up to the 
critical number of particles 
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where the factor 0.775 comes from the coefficient 
of the solution of the Lane-Emden equation with poly- 
tropic index n = 3, and mpi = is the Planck mass, 
mo is the nucleon mass and ji the average number of elec- 
trons per nucleon. As the kinetic energy of implosion is 
carried away by radiation the star settles down to a con- 
figuration of mass 



M = N clit m - U, 



(119) 



where the gravitational binding energy U can be as high 
as 5.72 x lO^Ncri^Q- 

Similarly Gamov B32I1 has shown that a gravitational 
collapse process to still higher densities can be stopped 
by the Fermi pressure of the neutrons (neutron star) and 
Oppenheimer ll33ll has shown that, if the effects of strong 
interactions are neglected, a configuration of equilibrium 
exists also in this case all the way up to a critical number 
of particles 

AU = 0.398^|, (120) 



where the factor 0.398 comes now from the integra- 
tion of the Tolman-Oppenheimer-Volkoff equation (see, 
e.g., Harrison et al. (1965) B3410 . If the kinetic energy 
of implosion is again carried away by radiation of pho- 
tons or neutrinos and antineutrinos the final configura- 
tion is characterized by the formula dl 19t with U < 
2.48 x 10 _2 A^ C iit'«o- These considerations and the exis- 
tence of such large values of the gravitational binding 
energy have been at the heart of the explanation of astro- 
physical phenomena such as red-giant stars and super- 
novae: the corresponding measurements of the masses of 
neutron stars and white dwarfs have been carried out with 
unprecedented accuracy in binary systems j35ll . 

From a theoretical physics point of view it is still an 
open question how far such a sequence can go: using 
causality nonviolating interactions, can one find a se- 
quence of braking and energy extraction processes by 
which the density and the gravitational binding energy 
can increase indefinitely and the mass-energy of the col- 
lapsed object be reduced at will? This question can also 
be formulated in the mass-formula language 10] (see also 
Ref. 14611 ): given a collapsing core of nucleons with a 
given rest mass-energy Mo, what is the minimum irre- 
ducible mass of the black hole which is formed? 

Following the previous two sections, consider a spher- 
ical shell of rest mass Mq collapsing in a flat space-time. 
In the neutral case the irreducible mass of the final black 
hole satisfies Eq. 11041 The minimum irreducible mass 
A/'™ 11 ' is obtained when the kinetic energy at the hori- 
zon T + is 0, that is when the entire kinetic energy T + has 
been extracted. We then obtain, form Eq. 11041 the simple 
result 

M (min) = Jfc (121) 



We conclude that in the gravitational collapse of a spher- 
ical shell of rest mass Mq at rest at infinity (initial energy 
M; = Mq), an energy up to 50% of Mqc 2 can in principle 
be extracted, by braking processes of the kinetic energy. 
In this limiting case the shell crosses the horizon with 
T + = 0. The limit ^ in the extractable kinetic energy 
can further increase if the collapsing shell is endowed 
with kinetic energy at infinity, since all that kinetic en- 
ergy is in principle extractable. 

We have represented in Fig. [10] the world lines of 
spherical shells of the same rest mass Mq, starting their 
gravitational collapse at rest at selected radii rg. These 
initial conditions can be implemented by performing 
suitable braking of the collapsing shell and concur- 
rent kinetic energy extraction processes at progressively 
smaller radii (see also Fig.fTTTi. The reason for the exis- 
tence of the minimum (1121b in the black hole mass is the 
"self closure" occurring by the formation of a horizon in 
the initial configuration (thick line in Fig.fTOb. 

Is the limit M; 1T — > ^ actually attainable without vio- 
lating causality? Let us consider a collapsing shell with 
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FIGURE 10. Collapse curves for neutral shells with rest 
mass Mq starting at rest at selected radii R* computed by 
using the exact solutions given in Ref. 1 47] . A different value 
of M m (and therefore of r + ) corresponds to each curve. The 
time parameter is the Schwarzschild time coordinate / and the 
asymptotic behaviour at the respective horizons is evident. The 
limiting configuration M- m = (solid line) corresponds to the 
case in which the shell is trapped, at the very beginning of its 
motion, by the formation of the horizon. 

charge Q. If M > Q an EMBH is formed. As pointed 
out in the previous section the irreducible mass of the fi- 
nal EMBH does not depend on the charge Q. Therefore 
Eqs. ( 1104b and ( 1121b still hold in the charged case. In 
Fig.QT|we consider the special case in which the shell is 
initially at rest at infinity, i.e. has initial energy M; = Mq, 
for three different values of the charge Q. We plot the 
initial energy M,, the energy of the system when all the 
kinetic energy of implosion has been extracted as well 
as the sum of the rest mass energy and the gravitational 

binding energy — ^ of the system (here ro is the radius 
of the shell). In the extreme case Q = Mq, the shell is 
in equilibrium at all radii (see Ref. [47]) and the kinetic 
energy is identically zero. In all three cases, the sum of 
the extractable kinetic energy T and the electromagnetic 

energy ^ reaches 50% of the rest mass energy at the 
horizon, according to Eq. ( 1121l i. 

What is the role of the electromagnetic field here? If 
we consider the case of a charged shell with Q ~ Mq, 
the electromagnetic repulsion implements the braking 
process and the extractable energy is entirely stored in 



the electromagnetic field surrounding the EMBH (see 
Ref. l46ll ). In the previous section we have outlined two 
different processes of electromagnetic energy extraction. 
We emphasize here that the extraction of 50% of the 
mass-energy of an EMBH is not specifically linked to 
the electromagnetic field but depends on three factors: 
a) the increase of the gravitational energy during the 
collapse, b) the formation of a horizon, c) the reduction 
of the kinetic energy of implosion. Such conditions are 
naturally met during the formation of an extreme EMBH 
but are more general and can indeed occur in a variety 
of different situations, e.g. during the formation of a 
Schwarzschild black hole by a suitable extraction of the 
kinetic energy of implosion (see Fig.[lO]and Fig.fTTI). 

Now consider a test particle of mass m in the gravi- 
tational field of an already formed Schwarzschild black 
hole of mass M and go through such a sequence of brak- 
ing and energy extraction processes. Kaplan [39] found 
for the energy E of the particle as a function of the radius 
r 

E = m ^\-m. (122) 

It would appear from this formula that the entire energy 
of a particle could be extracted in the limit r — > 2M. 
Such 100% efficiency of energy extraction has often been 
quoted as evidence for incompatibility between General 
Relativity and the second principle of Thermodynamics 
(see Ref. [40] and references therein). J. Bekenstein and 
S. Hawking have gone as far as to consider General Rel- 
ativity not to be a complete theory and to conclude that 
in order to avoid inconsistencies with thermodynamics, 
the theory should be implemented through a quantum 
description ll40ll42ll . Einstein himself often expressed the 
opposite point of view (see, e.g., Ref. 14111 and references 
therein). 

The analytic treatment presented can clarify this fun- 
damental issue. It allows to express the energy increase 
£ of a black hole of mass Mi through the accretion of a 
shell of mass Mq starting its motion at rest at a radius ro 
in the following formula which generalizes Eq. (1 1 22b : 

E = M-M X = - J| +M 0v /l - ^, (123) 

where M = M\ + E is clearly the mass-energy of the 
final black hole. This formula differs from the Kaplan 
formula ( 1122b in three respects: a) it takes into account 
the increase of the horizon area due to the accretion 
of the shell; b) it shows the role of the gravitational 
self energy of the imploding shell; c) it expresses the 
combined effects of a) and b) in an exact closed formula. 

The minimum value E m \ n of E is attained for the 
minimum value of the radius ro = 2M: the horizon of 
the final black hole. This corresponds to the maximum 




FIGURE 11. Energetics of a shell such that M\ = Mo, for selected values of the charge. In the first diagram Q = 0; the dashed 
line represents the total energy for a gravitational collapse without any braking process as a function of the radius R of the shell; 
the solid, stepwise line represents a collapse with suitable braking of the kinetc energy of implosion at selected radii; the dotted 
line represents the rest mass energy plus the gravitational binding energy. In the second and third diagram Q/Mq = 0.7, Q/Mq = 1 
respectively; the dashed and the dotted lines have the same meaning as above; the solid lines represent the total energy minus the 
kinetic energy. The region between the solid line and the dotted line corresponds to the stored electromagnetic energy. The region 
between the dashed line and the solid line corresponds to the kinetic energy of collapse. In all the cases the sum of the kinetic energy 
and the electromagnetic energy at the horizon is 50% of Mq. Both the electromagnetic and the kinetic energy are extractable. It is 
most remarkable that the same underlying process occurs in the three cases: the role of the electromagnetic interaction is twofold: 
a) to reduce the kinetic energy of implosion by the Coulomb repulsion of the shell; b) to store such an energy in the region around 
the EMBH. The stored electromagnetic energy is extractable as shown in Ref. 1491 



efficiency of the energy extraction. We have 
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or solving the quadratic equation and choosing the posi- 
tive solution for physical reasons 

E j ^ = 1(Jm 2 +M 2 -m\ (125) 

The corresponding efficiency of energy extraction is 

_ M -E mm _, 1 Mi I L , Ml A 
TJmax- Mq -1 2M lY i + ]gf I j U^>) 

which is strictly smaller than 100% for any given Mo ^ 
0. It is interesting that this analytic formula, in the limit 
M\ -c Mo, properly reproduces the result of equation 
dl211 >. corresponding to an efficiency of 50%. In the 
opposite limit M\ ^> Mo we have 

Tj max ~l-i^. (127) 

Only for Mo — > 0, Eq. ( 11261 ) corresponds to an efficiency 
of 100% and correctly represents the limiting reversible 
transformations. It seems that the difficulties of reconcil- 
ing General Relativity and Thermodynamics are ascrib- 
able not to an incompleteness of General Relativity but 
to the use of the Kaplan formula in a regime in which it 
is not valid. 



DYADOSPHERE FORMED IN 
GRAVITATIONAL COLLAPSES 

Electric field amplified in gravitational collapse. 

Initiating with supercritical electric fields on the core 
surface, we study pair production together with gravita- 
tional collapse. We use the exact solution of Einstein- 
Maxwell equations describing the gravitational collapse 
of a thin charged shell. Recall that the region of space- 
time external to the core is Reissner-Nordstrom with line 
element 

ds 2 = -a 2 dt 2 + a~ 2 dr 2 + r 2 dQ. 2 (128) 

in Schwarzschild like coordinate (t , r, 6 , 0) , where a 2 = 
1 — 2M/r + Q 2 /r 2 ; M is the total energy of the core as 
measured at infinity and Q is its total charge. Let us label 
with ro and to the radial and time-like coordinate of the 
core surface, and the equation of motion of the core is 

lima: 



Mo being the rest mass of the shell. The analytical solu- 
tions of Eq. ( 1129b were found to = to (ro) , and the core 
collapse speed V*(ro) as a function of ro is plotted in 
Fig. [12] where we indicate VI = V* \r =>- ds as the velocity 
of the core at the Dyadosphere radius r<j s . 

We now turn to the pair creation and plasma oscilla- 
tion taking place in the classical electric and gravitational 
fields during the gravitational collapse of a charged over- 
critical stellar core. As already show in Fig. [6] (i) the 
electric field oscillates with lower and lower amplitude 
around 0; (ii) electrons and positrons oscillates back and 
forth in the radial direction with ultra relativistic veloc- 
ity, as result the oscillating charges are confined in a thin 
shell whose radial dimension is given by the elongation 
Al = \l — /q I of the oscillations, where Iq is the radial co- 
ordinate of the center of oscillation and 

Al=[ A '^dt, (130) 
Jo pe 

where 7t e \\/p e = v is the radial mean velocity of charges. 
In Fig. [13] we plot the elongation Al as a function of 
time and electron mean velocity v as a function of the 
elongation during the first half period At of oscillation. 
This shows precisely the characteristic time At and size 
Al of charge confinement due to plasma oscillation. 

In the time At the charge oscillations prevent a macro- 
scopic current from flowing through the surface of the 
core. Namely in the time Af the core moves inwards of 

Ar* = V*Af>Al. (131) 

Since the plasma charges are confined within a region of 
thickness Al, due to Eq. ( 113U no charge "reaches" the 
surface of the core which can neutralize it and the initial 
charge of the core remains untouched. For example in the 
case M = 2OM ,^ = 0.1, and r = \r^, we have 

AZ<30Ac, Ar~10 3 T C , V* - 0.3c, (132) 

and Ar* ^> Al. We conclude that the core is not dis- 
charged or, in other words, the electric charge of the core 
is stable against vacuum polarization and electric field 
E = <2Ao is amplified during the gravitational collapse. 
As a consequence, an enormous amount (N ~ Q.r& s /eXc) 
(TBI as claimed in Ref. [4]) of pairs is left behind the col- 
lapsing core and Dyadosphere [4] is formed. 
Plasma expansion during gravitational collapse. 

The e + e~ pairs generated by the vacuum polariza- 
tion process around the core are entangled in the elec- 
tromagnetic field lHH], and thermalize in an electron- 
positron-photon plasma on a time scale ~ 10 4 Tc ifTlll 
(see Fig. [6]). As soon as the thermalization has occurred, 
the hydrodynamic expansion of this electrically neutral 
plasma starts lfl5l l45ll . While the temporal evolution of 
the e + e~y plasma takes place, the gravitationally col- 
lapsing core moves inwards, giving rise to a further am- 
plified supercritical field, which in turn generates a larger 
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FIGURE 12. Collapse velocity of a charged stellar core of mass Mq = 20M® as measured by static observers as a function of the 
radial coordinate of the core surface. Dyadosphere radii for different charge to mass ratios (£, = 10~ 3 , 1CT 2 , 10 _1 ) are indicated in 
the plot together with the corresponding velocity. 
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FIGURE 13. In left figure: Electrons elongation as function of time in the case r = The oscillations are damped in a time 

of the order of 10 3 — 10 4 Tc. The right figure: Electrons mean velocity as a function of the elongation during the first half oscillation. 
The plot summarize the oscillatory behaviour: as the electrons move, the mean velocity grows up from to the speed of light and 
then falls down at again. 



amount of e + e~ pairs leading to a yet higher temperature 
in the newly formed e + e~ y plasma. We report progress 
in this theoretically challenging process which is marked 
by distinctive and precise quantum and general relativis- 
tic effects. As presented in Ref. lfl3ll : we do not consider 
an already formed EMBH, but we follow the dynamical 
phase of the formation of Dyadosphere and of the asymp- 
totic approach to the horizon by examining the time vary- 
ing process at the surface of the gravitationally collaps- 
ing core. 

It is worthy to remark that the time-scale of hydrody- 
namic evolution (f ~ 0.1s) is, in any case, much larger 
than both the time scale needed for "all pairs to be 
created" (~ 10 3 Tc), and the thermalization time-scale 
(~ 10 4 Tc, see Fig. [6]) and therefore it is consistent to 
consider pair production, plus thermalization, and hydro- 
dynamic expansion as separate regimes of the system. 
We assume the initial condition that the Dyadosphere 
starts to be formed at the instant of gravitational col- 
lapse f(j s = to (r<j s ) = 0, and r<j s = R c the radius of mas- 
sive nuclear core. Having formulated the core collapse in 
General Relativity in Eq. ( 11291 ). we discretize the gravi- 
tational collapse of a spherically symmetric core by con- 
sidering a set of events (N— events) along the world line 
of a point of fixed angular position on the collapsing 
core surface. Between each of these events we consider 
a spherical shell of plasma of constant coordinate thick- 
ness Ar so that: 

1 . Ar is assumed to be a constant which is small with 
respect to the core radius; 

2. Ar is assumed to be large with respect to the mean 
free path of the particles so that the statistical de- 
scription of the e + e~ y plasma can be used; 

3. There is no overlap among the slabs and their union 
describes the entirety of the process. 

We check that the final results are independent of the 
special value of the chosen Ar and N. 

In each slab the processes of e + e~-pair production, 
oscillation with electric field and thermalization with 
photons are considered. While the average of the electric 
field S over one oscillation is 0, the average of S 2 is of 
the order of S 2 , therefore the energy density in the pairs 
and photons, as a function of ro, is given by 



£o(ro) = ^K 2 (ro)-<? c 2 ]=§ 



ro 



(133) 



For the number densities of e + e~ pairs and photons at 
thermal equilibrium we have n e + e - ~ tiy, correspond- 
ingly the equilibrium temperature To, which is clearly a 
function of ro and is different for each slab, is such that 



with e and n given by Fermi (Bose) integrals (with zero 
chemical potential): 
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From the conditions set by Eqs. ( 1134b . ( 1135b . ( 1136b . we 
can now turn to the dynamical evolution of the e + e~y 
plasma in each slab. We use the covariant conservation of 
energy momentum and the rate equation for the number 
of pairs in the Reissner-Nordstrom geometry external to 
the core: 



V a T ab = 0, 



(137) 



V a {n e+e -u a ) = ov [n 2 e+e - (T) - « 2 +e _] , (138) 



where T a0 = {£ + p)u"u D + pg ab is the energy- 
momentum tensor of the plasma with proper energy 
density e and proper pressure p, u" is the fluid 
4— velocity, n e + e - is the number of pairs, n e + e - (T) 
is the equilibrium number of pairs and ffv is the mean of 
the product of the e + e~ annihilation cross-section and 
the thermal velocity of pairs. In each slab the plasma 
remains at thermal equilibrium in the initial phase of the 
expansion and the right hand side of the rate Eq .J[138b 
is effectively 0, see Fig. 24 (second panel) of 14911 for 
details. 

If we denote by % a the static Killing vector field nor- 
malized at unity at spacial infinity and by {E r } ( the fam- 
ily of space-like hypersurfaces orthogonal to it being 
the Killing time) in the Reissner-Nordstrom geometry, 
from Eqs. (1138b . the following integral conservation laws 
can be derived 



E, a T ab dY. h =E, / n e+e -u b dL h =N e+e -, (139) 



where dT./, = a _2 ^r 2 sin QdrdQdtj) is the vector surface 
element, E the total energy and N e + e - the total number of 
pairs which remain constant in each slab. We then have 



[(e+p)y 2 -p]r 1 = <£, n e+e -ya- l r 2 =^ e 



where £ and OT^+g- are constants and 



y=crV§ fl = i-cr 4 (£) 
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e(Zb) = e r (T ) + e e+ (T )+e e - (T ) = e , (134) 



is the Lorentz 7 factor of the slab as measured by static 
observers. We can rewrite Eqs. (1139b for each slab as 



(£) =« 4 / ro , 



(142) 

/'o-l-(^) 2 (^) 2 fe) 4 (144) 



where pedex refers to quantities evaluated at se- 
lected initial times to > 0, having assumed r(?o) = ro, 
dr/dt\ 1=to =0,T(t ) = T . 

Eq. (1142b is only meaningful when f rg (r) > 0. From 
the structural analysis of such equation it is clearly iden- 
tifiable a critical radius ro such that: 

• for any slab initially located at ro > R we have 
fr (r) > for any value of r > ro and /,-„ (r) < for 
r < ro; therefore a slab initially located at a radial 
coordinate ro > R moves outwards, 

• for any slab initially located at ro < R we have 
fr {r) > for any value of r + < r < ro and /,-„ (r) < 
for r > ro; therefore a slab initially located at 
a radial coordinate ro < R moves inwards and is 
trapped by the gravitational field of the collapsing 
core. 

We define the surface r — R, the Dyadosphere trap- 
ping surface (DTS). The radius R of DTS is generally 

evaluated by the condition = 0. R is so close 

to the horizon value r+ that the initial temperature Tq 
satisfies kTo ^> m e c 2 and we can obtain for R an ana- 
lytical expression. Namely the ultra relativistic approx- 
imation of all Fermi integrals, Eqs. ( 1135b and ( 1136b , 
is justified and we have n e + e - (T) °< T 3 and therefore 
f ro ^l- (T/T0) 6 (Oo/a) 2 (r/r ) 4 (r < R). The denning 
equation of R, together with (1144b . then gives 



R = 2M 1 + (1 



■3Q 2 /4M 2 ) l/2 } >r 4 



(145) 



In the case of an EMBH with M = 2OM , Q = 0. 1M, we 
compute: 

• the fraction of energy trapped in DTS: 



E = / ae dT, ~ 0.53 / ae Q dI.; 

Jr + <r<R Jr + <r<r is 

(146) 

• the world-lines of slabs of plasma for selected ro in 
the interval (R, r^) (see left figure in Fig. [141 : 

• the world-lines of slabs of plasma for selected ro in 
the interval {r+,R) (see Fig.fTBTl. 

At time f = ?o (R) when the DTS is formed, the plasma 
extends over a region of space which is almost one order 
of magnitude larger than the Dyadosphere and which 



we define as the effective Dyadosphere. The values of 
the Lorentz 7 factor, the temperature and e + e~ number 
density in the effective Dyadosphere are given in the 
right figure in Fig. [14] 



HYDRODYNAMIC EXPANSION AFTER 
GRAVITATIONAL COLLAPSE. 

Plasma expansion after gravitational collapse. 

After gravitational collapse, the adiabatically hydro- 
dynamic expansion of e + e~ 7 plasma continues, obeying 
hydrodynamic and rate equations ( 11371138b for conser- 
vations of energy-momentum, entropy and particle num- 
bers, until it becomes optically thin: 



dr(n e ±)a T ^0{\), 



(147) 



where Oj — 0.665 x 10 cm is the Thomson cross- 
section and the integration is over the radial interval of 
the plasma in the comoving frame. At this point the en- 
ergy is virtually entirely in the form of free-streaming 
photons. The calculations were independently performed 
by numerical simulation in Lawrence Livermore Na- 
tional Laboratory using the one dimensional (1-D) hy- 
drodynamic code, and analytical approach in ICRA, Uni- 
versity of Rome lfl5l l45ll using approximate hydrody- 
namic and rate equations d 1 371 1 38b neglecting gravita- 
tional effects, 
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where the thermal index r=l+p/e,^ / is the volume of 
a single slab, N e + e - = Jn e + e - is the pair number density 
as measured in the laboratory frame by an observer at 
rest with the black hole, and N e + e - (T) is the equilibrium 
laboratory pair number density. In Fig. [16] we show the 
Lorentz gamma factor as a function of radius and the 
decoupling gamma factor at the transparency point dl47b 
as a function of EMBH masses. We find the expansion 
pattern of a shell of constant coordinate thickness with 
the astrophysically unprecedented large Lorentz factors, 
named by the Pair-ElectroMagnetic Pulse (PEM-pulse). 
Pair- annihilation and reheating. 

In Fig. [17] we plot the number-densities of pairs n e + e - 
given by the rate equation ( 1 149b and n e + e - (T) computed 
from a Fermi-integral with zero chemical potential (1136b 
with the temperature T determined by the equilibrium 
condition (1134b . It clearly indicates that the pairs e 
fall out of equilibrium as the temperature drops below 
the threshold of e ± -pair annihilation. As a consequence 
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FIGURE 14. In left figure: World line of the collapsing charged core (dashed line) as derived from Eq. d 129b : world lines of 
slabs of plasma for selected radii ro in the interval (R, r^ s ). At time t the expanding plasma extends over a region which is almost 
one order of magnitude larger than the Dyadosphere. The small rectangle in the right bottom is enlarged in Fig. Q3] The right 
figure: Physical parameters in the effective Dyadosphere: Lorentz / factor, proper temperature and proper e + e~ number density as 
functions at time t. 



of pair e -annihilation, the crossover of this reheating 
process is also shown in Fig. [17] This can be understood 
as follows. From the conservation of entropy it follows 
that asymptotically we have 



{VT 



3 )r<mt 2 _ 11 



(vr 3 ) 



T>mc 



(150) 



where the comoving volume V = ~f /y. The same consid- 
erations when repeated for the conservation of the total 
energy eyV — ey 2 f / following from Eq. ( 1148b then lead 
to 

(VT 4 y) T<mc 2 = n 

4 



(151) 



(VT 4 y) T>mc 2 

The ratio of these last two relations then gives asymptot- 
ically 

To = (Ty) T>mc 2 = (Ty) T<mc 2, (152) 

where T a is the initial average temperature of the Dyado- 
sphere at rest, given in Fig. [2] Eq. ( 11481 ) also explains the 
approximate constancy of Ty shown in Fig.fTSI 
Photon temperature at transparency. 

We are interested in the observed spectrum at the time 
of decoupling. To calculate the spectrum, we assume that 
(i) the plasma fluid of coupled e + e~ pairs and photons 
undergoes adiabatic expansion and does not emit radia- 
tive energy before they decouple f50ll : (ii) the e + <?~ pairs 



and photons are in equilibrium at the same temperature T 
when they decouple. Thus the photons that are described 
by a Planck distribution in an emitter's rest-frame with 
temperature T' will appear Planckian to a moving ob- 
server, but with boosted temperature T 



(6,v,T') 



r 



exp( 



1 



y(l-^cos0)' 
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where ~ cos 9 is the component of plasma velocity di- 
rected toward the observer. We integrate over angle with 
respect to the observer, and get the observed number 
spectrum N e , per photon energy e, per steradian, of a 
relativistically expanding spherical shell with radius R, 
thickness dR in cm, velocity v, Lorentz factor y and fluid- 
frame temperature T' to be (in photons/eV/4?r) 



N e (v,T',R)= dV 



eT' 

(5.23 x \0 n )AnR 2 dR x 

vy 



x log 



l-exp[-ye(l + j;)/r / ] 
l-«p[-yB(l-2)/r'] 



(154) 



which has a maximum at e mux = 1.39/T' eV for y^> 1. 
We may then sum this spectrum over all shells of our 
PEM-pulse to get the total spectrum. Since we had as- 
sumed the photons are thermal in the comoving frame, 
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FIGURE 15. Enlargement of the small rectangle in the right bottom of left figure in Fig.QJ] World-lines of slabs of plasma for 
selected radii ro in the interval (r + R). 



our spectrum (1154b has now an high frequency expo- 
nential tail, and the spectrum appear as not thermal. 
Fig. [18] illustrates the extreme relativistic nature of the 
PEM pulse expansion: at decoupling, the local comoving 
plasma temperature is < 1 keV, but is boosted by yfi na i 
to ~ 1 MeV and its relation to the initial energy of the 
PEM pulse in the Dyadosphere is presented. In Ref. i5lll 
we obtain the observed light curve by decomposing the 
spherical PEM-pulse into concentric shells, considering 
the light curve in the relative arrival time t a of the first 
light from each shell, then summing over contributions 
from all shells. 



PREDICATIONS IN CONNECTION 
WITH SHORT GAMMA RAY BURSTS 

The formation of Dyadosphere during gravitational 
collapse and its hydrodynamic expansion after gravita- 
tional collapse continuously connect when gravitational 
collapse ends. Eqs. d 14-81 149b must be integrated 
and matched with the solution found by integrating 
Eqs. d 1 371 1 38t at the transition between the two regimes. 
In this general framework we analyze in Refs. |[52l l53ll . 
the gravitational formation and then hydrodynamic 
evolution of the Dyadosphere. We recall that separatrix 
was found in the motion of the plasma at a critical 
radius R (11451 ): the plasma created at radii smaller than 
R is trapped by the gravitational field of the collapsing 
core and implodes toward the black hole, while the one 
created at radii larger than R expands outward. The 
plasma (r > R) is divided into N slabs with thickness 
Ar, to describe the adiabatic hydrodynamic expansion 
of the optically thick plasma slabs (the PEM-pulse) all 
the way to the point where the transparency condition 



is reached. Eqs. ( fT37b . (HD and ( 11481149b have been 
integrated and results are presented in Fig. [19] The 
integration stops when each slab of plasma reaches the 
optical transparency condition given by 




where the integral extends over the radial thickness Ar of 
the slab. The overall independence of the result of the dy- 
namics on the number of the slabs adopted in the dis- 
cretization process or analogously on the value of Ar has 
also been checked. We have repeated the integration for 
N = 10, N = 100 reaching the same result to extremely 
good accuracy. The results in Fig. [T9l correspond to the 
case N =10. The evolution of each slab occurs without 
any collision or interaction with the other slabs; see the 
left diagram in Fig. [19] The outer layers are colder and 
less dense than the inner ones and therefore reach trans- 
parency earlier; see the right diagram in Fig. [19] 

As the plasma becomes transparent, gamma ray pho- 
tons are emitted, and the luminosity and spectrum are 
calculated, analogously to Eqs. (11531154b . In Fig. [20] 
where we plot both the theoretically predicted luminos- 
ity L and the spectral hardness r o b s of the signal reach- 
ing a far-away observer as functions of the arrival time 
t a . Since three quantities depend in an essential way on 
the cosmological redshift factor z, see Refs. IBH l54Tl . 
we have adopted a cosmological redshift z = 1 for this 
figure. The energy of the observed photon is kT^ = 
kjT' j (1 +z), where k is the Boltzmann constant, T' is 
the temperature in the comoving frame of the pulse and 
7 is the Lorentz factor of the plasma at the transparency 
time. The initial zero of time is chosen as the time when 
the first photon is observed, then the arrival time t a of a 
photon at the detector in spherical coordinates centered 
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FIGURE 16. In left figure: Lorentz gamma factor y as a function of radius. Three models for the expansion pattern of the PEM- 
pulse are compared with the results of the one dimensional hydrodynamic code for a IOOOMq black hole with charge to mass ratio 
% =0.1. The 1-D code has an expansion pattern that strongly resembles that of a shell with constant coordinate thickness. The right 
figure: In the expansion model of a shell with constant coordinate thickness, the decoupling gamma-factor yg na i as a function of 
EMBH masses is plotted with charge to mass ratio t, =0.1. 
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FIGURE 17. In left figure: The number-densities (A^i/cm -3 ) of pairs n e ± (solid line) as obtained from the rate equation d 149b , 
and n e ± (T) (dashed line) as computed by Fermi-integrals with zero chemical potential d 1 36b . provided the temperature T determined 
by the equilibrium condition d 1 34b . are plotted for a IOOOMq black hole with charge to mass ratio S, =0.1. For T < m e c 2 , the two 
curves strongly diverge. The right figure: The reheating ratio T^V /T^V a defined by Eq. ( 1 1 50b is plotted as a function of radius for a 
1000M ;, black hole with charge to mass ratio £, =0.1. The rate equation d 1 49b naturally leads to the value -j- after e + e~ -annihilation 
has occurred. 
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FIGURE 18. In left figure: Temperature of the plasma is shown over a range of masses. 7j n i t is the average initial temperature of 
the plasma deposited around the EMBH, T \, s = 1.39yfi na ir' is the observed peak temperature of the plasma at decoupling while T' 
is the comoving temperature at decoupling. Notice that yT ~ 7j n i t as expected from Eq. ( 1152b . The right figure: The peak of the 
observed number spectrum as a function of the EMBH mass is plotted with charge to mass ratio | =0.1. 
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FIGURE 19. Expansion of the plasma created around an overcritical collapsing stellar core with M = WM,., and Q = 0. \\J~GM. 
Left diagram: world lines of the plasma. Right diagram: Lorentz y factor as a function of the radial coordinate r, showing the 
PEM-pulse reaches ultra relativistic regimes with Lorentz factor y ~ 10 2 -10 4 . 




FIGURE 20. Predicted observed luminosity and observed 
spectral hardness of the electromagnetic signal from the grav- 
itational collapse of a collapsing core with M = IOMq, Q = 
0.1 yGM at z = 1 as functions of the arrival time t a . 



on the black hole is given by I5lll54ll : 

ta = (l+z) 



f + £0 _^ cos0 



(156) 



where (t,r(t) ,9,<j>) labels the laboratory emission event 
along the world line of the emitting slab and ro is the 
initial position of the slab. The projection of the plot in 
Fig. [20] onto the t u -L plane gives the total luminosity as 
the sum of the partial luminosities of the single slabs. 
The sudden decrease of the intensity at the time f„ = 
0.040466 s corresponds to the creation of the separatrix 
introduced in Ref. l44ll . We find that the duration of 
the electromagnetic signal emitted by the relativistically 
expanding pulse is given in arrival time by 



5 x 10~ 2 s. 



(157) 



The projection of the plot in Fig. [20] onto the kT \, s , t a 
plane describes the temporal evolution of the spectral 
hardness. We observe a precise soft-to-hard evolution of 
the spectrum of the gamma ray signal from ~ 10 2 KeV 
monotonically increasing to ~ 1 MeV. The correspond- 
ing spectra are presented in Figs.[2T] which are consistent 
with observations of short GRBs 15311 . The three quan- 
tities L,kT \, s and t a are clearly functions of the charge 



Q and the mass M of the collapsing core. We present 
in Fig. [22] the arrival time interval for M ranging from 
M ~ 1OM to 1O 3 M , keeping Q = 0.1 VGM. The arrival 
time interval is very sensitive to the mass of the black 
hole: 

A* a ~ lO^-lO^s. (158) 



Similarly the spectral hardness of the signal is sensitive 
to the ratio Q/VGM f52ll . Moreover the duration, the 
spectral hardness and luminosity are all sensitive to the 
cosmological redshift z (see Ref. lf52ln . 

Ghirlanda et al. have given evidence for the existence 
of an exponential cut off at high energies in the spectra of 
short GRBs. We are currently comparing and contrasting 
these observational results with the predicted cut off in 
Fig. [20]which results from the existence of the separatrix 
introduced in ill 311 . The observational confirmation of the 
results presented in Fig. [20]would lead for the first time 
to the identification of a process of gravitational collapse 
and its general relativistic self-closure as seen from an 
asymptotic observer. 

The characteristic spectra, time variabilities and lumi- 
nosities of the electromagnetic signals from collapsing 
overcritical stellar cores, here derived from first princi- 
ples, a gree s very closely with the observations of short 
GRBs [55]. New space missions must be planned, with 
temporal resolution down to fractions of /is and higher 
collecting area and spectral resolution than at present, in 
order to verify the detailed agreement between our model 
and the observations. It is now clear that if these theo- 
retical predictions will be confirmed, we would have a 
very powerful tool for cosmological observations: the in- 
dependent information about luminosity, time-scale and 
spectrum can uniquely determine the mass, the electro- 
magnetic structure and the distance from the observer of 
the collapsing core, see, e.g., Fig. [22] and Ref. l52l l53ll . 
In that case short-bursts may become the best example of 
standard candles in cosmology. 
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